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ABSTRACT 

The  brightest  cluster  radio  halo  known  resides  in  the  Coma  cluster  of  galaxies.  The  relativistic  electrons  producing 
this  diffuse  synchrotron  emission  should  also  produce  inverse  Compton  emission  that  becomes  competitive  with 
thermal  emission  from  the  intracluster  medium  (ICM)  at  hard  X-ray  energies.  Thus  far,  claimed  detections  of  this 
emission  in  Coma  are  controversial.  We  present  a  Suzaku  HXD-PIN  observation  of  the  Coma  cluster  in  order  to 
nail  down  its  nonthermal  hard  X-ray  content.  The  contribution  of  thermal  emission  to  the  HXD-PIN  spectrum 
is  constrained  by  simultaneously  fitting  thermal  and  nonthermal  models  to  it  and  a  spatially  equivalent  spectrum 
derived  from  an  XMM-Newton  mosaic  of  the  Coma  field.  We  fail  to  find  statistically  significant  evidence  for 
nonthermal  emission  in  the  spectra  which  are  better  described  by  only  a  single-  or  multitemperature  model  for 
the  ICM.  Including  systematic  uncertainties,  we  derive  a  90%  upper  limit  on  the  flux  of  nonthermal  emission  of 
6.0  x  10-12  erg  s-1  cm-2  (20-80  keV,  for  T  =  2.0),  which  implies  a  lower  limit  on  the  cluster-averaged  magnetic 
field  of  B  >  0.15  /xG.  Our  flux  upper  limit  is  2.5  times  lower  than  the  detected  nonthermal  flux  from  RXTE  and 
BeppoSAX.  However,  if  the  nonthermal  hard  X-ray  emission  in  Coma  is  more  spatially  extended  than  the  observed 
radio  halo,  the  Suzaku  HXD-PIN  may  miss  some  fraction  of  the  emission.  A  detailed  investigation  indicates  that 
~50%-67%  of  the  emission  might  go  undetected,  which  could  make  our  limit  consistent  with  that  of  Rephaeli  & 

Gruber  and  Fusco-Femiano  et  al.  The  thermal  interpretation  of  the  hard  Coma  spectrum  is  consistent  with  recent 
analyses  of  INTEGRAL  and  Swift  data. 
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1.  INTRODUCTION 

In  the  hierarchical  scenario  of  cosmic  structure  formation, 
clusters  of  galaxies  form  at  late  times  through  mergers  between 
subclusters  and  through  the  accretion  of  galaxies  and  galaxy 
groups.  The  distribution  of  their  massive  halos  in  space  and 
time  depend  sensitively  on  the  underlying  cosmology,  and 
much  effort  has  been  made  to  connect  observable  properties 
of  the  gas  to  the  total  cluster  mass  in  order  to  constrain 
cosmological  parameters  (e.g.,  Mantz  et  al.  2008,  and  references 
therein).  However,  merger  processes  are  known  to  significantly 
disrupt  the  thermal  gas  (e.g.,  Ricker  &  Sarazin  2001  Ritchie 
&  Thomas  2002),  typically  biasing  inferred  masses  and  the 
resultant  cosmological  parameter  estimates  (Randall  et  al.  2002; 
Wik  et  al.  2008).  Merger-induced  shocks  and  turbulence,  besides 
heating  the  gas,  are  thought  to  also  reaccelerate  relativistic 
particles  present  in  the  intracluster  medium  (ICM;  Sarazin 
1999;  Brunetti  &  Blasi  2005).  Nonthermal  electrons,  observed 
via  diffuse,  radio  synchrotron  emission,  have  been  detected  in 
over  50  clusters,  all  of  them  undergoing  mergers  (Buote  2001; 
Schuecker  et  al.  2001).  If  the  energy  in  a  relativistic  phase  of  the 
ICM  is  large  enough  to  add  pressure  support  to  the  thermal  gas, 
even  transiently,  the  ability  to  derive  masses  and  therefore  use 
clusters  as  cosmological  probes  may  be  compromised  (Skillman 
et  al.  2008).  An  assessment  of  the  relativistic  contribution  to  the 


energy  budget  of  clusters  is  necessary  to  fully  characterize  the 
state  of  the  ICM. 

Diffuse,  cluster- wide  synchrotron  radio  emission,  called  radio 
halos  or  relics  depending  on  their  morphology,  imply  that  both 
magnetic  fields  and  relativistic  electron  populations  are  present 
on  large  scales.  The  total  luminosity  of  a  synchrotron-emitting 
electron  is  given  by 

4  9 

Lr  =  - crTcy  €B ,  (1) 

where  aT  is  the  Thomson  cross  section,  c  is  the  speed  of  light, 
y  is  the  Lorentz  factor  of  the  electron,  and  eB  =  B2/%tt  is  the 
energy  density  of  the  magnetic  field.  For  many  such  electrons, 
the  value  of  LR  depends  both  on  the  number  of  electrons  and  on 
B  and  cannot  independently  determine  either.  However,  these 
same  electrons  will  up- scatter  cosmic  microwave  background 
(CMB)  photons  through  inverse  Compton  (IC)  interactions, 
which  have  a  luminosity  Lx  equivalent  in  form  to  Equation  (1) 
but  with  €B  replaced  by  the  energy  density  of  the  CMB.  Since 
both  luminosities  are  proportional  to  the  number  of  electrons, 
their  ratio  gives  the  volume- averaged  magnetic  field, 

lr  _  B2/8n 

LX  d  1  CMB 
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where  a  is  the  radiation  constant  and  Tcmb  is  the  temperature 
of  the  CMB.  The  IC  radiation  should  be  observable  at  hard  X- 
ray  energies  (Rephaeli  1977).  Thus  far,  IC  emission  has  only 
been  detected  at  low  significance  (Nevalainen  et  al.  2004)  or,  in 
one  case,  in  a  cluster  with  weaker  radio  emission  (Eckert  et  al. 
2008a;  but  see  also  Ajello  et  al.  2009  and  Fujita  et  al.  2008).  The 
measurement  of  an  IC  flux  from  a  synchrotron  source  directly 
leads  to  a  simultaneous  determination  of  the  average  value  of  B 
and  the  relativistic  electron  density  (Harris  &  Romanishin  1974; 
Sarazin  1988).  Therefore  searches  for  IC  emission  coincident 
with  a  radio  halo  or  relic  are  an  excellent  way  to  constrain  the 
contribution  of  relativistic  materials  in  clusters. 

The  first,  and  brightest,  radio  halo  was  discovered  by  Willson 
(1970)  in  the  Coma  cluster,  and  its  radio  properties  have  perhaps 
been  the  best  studied  (e.g.,  Giovannini  et  al.  1993  Deiss  et  al. 
1997;  Thierbach  et  al.  2003).  Coma  has  been  observed  by  all 
the  major  observatories  with  hard  X-ray  capabilities  (Rephaeli 
et  al.  1994;  Hughes  et  al.  1993;  Bazzano  et  al.  1990;  Henriksen 
&  Mushotzky  1986),  and  more  recently  nonthermal  detections 
have  been  claimed  by  Rephaeli  &  Gruber  (2002)  with  RXTE  and 
by  Fusco-Femiano  et  al.  (1999,  2004)  with  BeppoSAX ,  though 
the  latter  detection  is  controversial  (Rossetti  &  Molendi  2004; 
Fusco-Femiano  et  al.  2007).  Due  to  the  large  field  of  view  (FOV) 
of  these  nonimaging  instruments  and  the  simple  characterization 
of  the  thermal  gas,  the  source  of  this  emission  remains  uncertain. 
Even  more  recently,  long  (~  1  Ms)  observations  with  INTEGRAL 
have  imaged  extended  diffuse  hard  X-ray  emission  from  Coma, 
though  it  was  found  to  be  completely  consistent  with  thermal 
emission  (Renaud  et  al.  2006;  Eckert  et  al.  2007;  Lutovinov 
et  al.  2008). 

In  this  study,  we  present  a  Suzaku  HXD-PIN  observation  of 
the  Coma  cluster  in  an  effort  to  detect  nonthermal  emission  as¬ 
sociated  with  the  radio  halo  and  potentially  confirm  the  RXTE 
and  BeppoSAX  detections.  The  HXD-PIN  instrument  has  a  non¬ 
imaging  collimator  like  those  on-board  RXTE  and  BeppoSAX , 
but  with  a  FOV  about  a  quarter  as  large,  which  reduces  the  pos¬ 
sible  contamination  from  hard  point  sources  (Takahashi  et  al. 
2007).  Also,  the  Suzaku  particle  background  is  ~5  times  lower 
than  the  backgrounds  of  either  RXTE  or  BeppoSAX  (Mitsuda 
et  al.  2007).  In  order  to  clearly  distinguish  the  thermal  and  non¬ 
thermal  emission  visible  within  the  PIN,  the  hard  Suzaku  PIN 
spectrum  is  jointly  fit  with  a  spatially  equivalent  XMM-Newton 
EPIC-pn  spectrum.  The  XMM  spectrum,  at  lower  energies  and 
completely  dominated  by  thermal  emission,  allows  Coma’s 
thermal  and  nonthermal  properties  to  be  simultaneously  de¬ 
termined.  The  XMM  and  HXD-PIN  observations  are  reported  in 
Section  2  and  the  extraction  of  the  resulting  spectra  is  discussed 
in  Section  3.  Fits  to  the  joint  spectra  are  described  in  Section  4. 
In  Section  5,  we  discuss  the  implications  of  our  results  for  the 
nature  of  the  hard  X-ray  emission  from  the  Coma  cluster.  We 
assume  a  flat  cosmology  with  QM  =  0.23  and  Hq  =  72  km  s-1 
Mpc-1  and  a  luminosity  distance  to  Coma  of  98.4  Mpc.  Unless 
otherwise  stated,  all  uncertainties  are  given  at  the  90%  confi¬ 
dence  level. 

2.  OBSERVATIONS 

The  Suzaku  observation  was  undertaken  as  part  of  AO-1  from 
2006  May  31  through  June  4,  soon  after  16  of  the  64  PIN  diode 
bias  voltages  were  lowered  from  500  V  to  400  V,  but  before 
an  additional  16  diodes  were  similarly  lowered.  We  analyze 
Version  2  of  the  pre-processed  data  (PROCVER  2.0.6.13), 
which  allows  for  the  diode  bias  drop,  with  HEAsoft  6.4.0  and 
XSPEC  12.4.0w.  For  the  HXD-PIN  instrument,  the  standard 


data  selection  criteria  are  applied  to  extract  the  source  spectrum, 
and  the  same  criteria  are  used  to  select  times  for  the  modeled 
non-X-ray  background  (NXB)  spectrum.  Specifically,  we  select 
observing  times  when  the  geomagnetic  cut-off  rigidity  (COR) 
is  above  a  critical  value  (COR  >  6  GV),  when  the  satellite  is 
not  within  the  South  Atlantic  Anomaly  (SAA_HXD  =  0)  or 
has  just  left  it  (T_SAA_HXD  >  500  s),  and  when  Suzaku  is 
pointed  above  and  at  least  5°  away  from  the  Earth’s  horizon 
(ELV  >  5°).  The  strength  of  the  NXB  is  known  to  be  roughly 
inversely  proportional  to  the  value  of  the  COR  and  to  be  elevated 
inside  the  SAA,  gradually  decaying  to  typical  levels  after  SAA 
passage.  These  criteria  ensure  that  the  low  NXB  of  the  HXD  is 
minimized  and  can  be  well  characterized,  which  is  necessary 
if  it  is  to  be  accurately  modeled.  After  event  selection,  the 
PIN  exposure  time  is  reduced  from  166.2  ks  to  156.1  ks  after 
dead-time  correction.  The  HXD-GSO  spectrum  is  found  to  be 
consistent  with  the  background,  so  we  do  not  consider  it  further 
here.  We  use  the  Suzaku  CCD  data  from  the  XIS0  chip  to  check 
the  cross-calibration  of  Suzaku  and  XMM-Newton.  Standard 
event  selection  was  applied  to  the  XIS0  data,  leading  to  an 
exposure  time  of  178.7  ks. 

The  mosaic  XMM-Newton  observations  of  the  Coma  cluster, 
including  14  separate  pointings,  were  done  as  part  of  an 
instrument  performance  verification  program,  a  complete  log 
of  which  is  presented  in  Schuecker  et  al.  (2004).  The  initial 
observations  were  undertaken  by  and  first  reported  in  Briel 
et  al.  (2001).  Seven  new  observations,  aimed  at  resolving 
the  temperature  structure  of  the  Coma  center,  have  also  been 
performed  (PI  P.  Schuecker).  However,  high  solar  activity 
during  the  exposures  resulted  in  a  high  detector  background 
above  2  keV,  making  these  observations  less  suitable  for  our 
purposes,  and  therefore  we  use  only  the  observations  reported  in 
Section  4  and  Table  2  of  Schuecker  et  al.  (2004).  We  choose  only 
to  include  the  EPIC-pn  data  from  XMM  in  our  analysis.  Because 
these  observations  were  made  early  in  the  mission,  they  cannot 
be  processed  with  the  standard  software,  though  the  EPIC-pn 
data  have  undergone  in-house  processing.  Also,  its  effective 
area  at  high  energies  is  higher  than  for  the  EPIC-MOS  detector, 
making  it  the  more  suitable  instrument.  The  benefit  of  including 
the  EPIC-MOS  data  is  unclear,  due  to  the  addition  of  cross¬ 
calibration  errors  and  given  the  already  high  signal  to  noise  of 
the  EPIC-pn  data. 

3.  EXTRACTION  AND  CONSTRUCTION  OF  SPECTRA 

To  produce  complementary  spectra  from  the  XMM  EPIC-pn 
and  Suzaku  HXD-PIN  data  that  can  be  simultaneously  fit,  the 
background  and  responses  of  both  instruments  must  be  carefully 
considered  to  minimize  systematic  uncertainties.  The  expected 
nonthermal  signal  is  near  the  limit  of  the  PIN  sensitivity,  and  a 
robust  characterization  of  this  emission  particularly  depends  on 
the  accuracy  of  the  PIN  background  and  XMM-Suzaku  cross¬ 
normalization. 

3.1.  HXD-PIN  Spectrum  and  Non-X-Ray  Background 

As  the  HXD  is  a  nonimaging  instrument,  we  simply  extract 
the  PIN  spectrum  from  the  selected  events  and  group  the  spectral 
bins  so  that  each  bin  contains  at  least  30  counts  to  ensure  that 
Gaussian  statistics  and  x2  fitting  are  valid.  The  response  matrix 
is  provided  in  the  Suzaku  CALDB8  for  Version  2  data  products, 
and  we  use  ae_hxd_pinhxnome2_20080129.rsp  for  all  source 
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Figure  1.  Suzaku  HXD-PIN  NXB  model  spectrum  (red  data  points)  compared 
with  the  Coma  cluster  data  (uncorrected  for  background).  Note  that  at  energies 
above  45  keV,  the  NXB  dominates  the  data  and  that  deviations  of  the 
data  above  the  NXB  are  confined  to  individual  channels  that  are  simply 
statistical  fluctuations  or  are  otherwise  imperfectly  characterized  by  the  NXB 
model. 

components  other  than  the  cosmic  X-ray  background  (CXB), 
for  which  ae_hxd_pinflate2_20080129.rsp  is  used. 

The  non-X-ray  background  for  a  PIN  observation  is  most 
accurately  obtained  from  a  model,  as  opposed  to  a  comparable 
blank  field  observation.  This  method  is  motivated  by  the  strong 
dependence  of  the  background  count  rate  and  spectral  shape 
on  the  value  of  the  geomagnetic  COR  and  the  time  since  the 
passage  of  the  satellite  through  the  South  Atlantic  Anomaly 
(SAA),  quantities  which  vary  and  have  a  unique  distribution  for 
every  observation.  The  model  matches  the  distribution  of  the 
COR  and  other  parameters  of  the  observed  data.  A  model  is 
also  required  because  there  is  no  concurrent  measurement  of 
the  NXB,  such  as  by  nodding  between  Coma  and  a  blank  field. 
We  use  the  so-called  “bgd_d”  model  for  Version  2  processed 
data,  which  makes  use  of  HXD-GSO  information  as  well  as 
the  COR  and  SAA  values.  This  NXB  model  is  shown  with  the 
PIN  data  spectrum,  uncorrected  for  the  background,  in  Figure  1 . 
While  the  shape  of  the  NXB  is  generally  well  reproduced,  the 
success  of  the  model  in  determining  its  overall  normalization  is 
±2.3%  from  15-40  keV  and  ±4%  from  40-70  keV  (Mizuno 
et  al.  2008).  These  estimates  of  the  systematic  error  in  the  NXB 
are  extrapolated  to  the  90%  confidence  interval  from  the  la 
values  derived  using  Earth-occulted  data  in  Mizuno  et  al.  (2008). 
We  adopt  these  values  (2.3%  from  12-40  keV  and  4%  from  40- 
70  keV)  as  our  estimate  of  the  90%  systematic  error  in  the  PIN 
NXB.  To  confirm  the  accuracy  of  the  model  background,  we 
extracted  events  for  both  the  data  and  model  for  times  when 
the  Earth  occults  the  PIN  FOV  (ELV  <  —5°,  all  other  selection 
criteria  unchanged);  the  resulting  spectra  are  shown  in  Figure  2. 
The  fractional  difference  between  the  model  NXB  count  rate  and 
the  Earth-occulted  data  is  (0.2  ±  1.0)%  over  the  energies  12- 
40keV  and(— 2.0 ±2.6)%  from  40-70  keV  (lcr  errors).  Over  the 
whole  range  considered,  the  fractional  difference  in  count  rates 
is  extremely  small:  (0.005  ±  0.9)%.  Because  of  the  excellent 
agreement,  we  do  not  adjust  the  level  of  the  background  as 
proposed  in  Ishida  et  al.  (2007a). 

The  estimate  of  the  systematic  error  adopted  here  comes  from 
an  analysis  of  Earth-occulted  data,  which  is  the  same  data  used 
to  generate  models  of  the  NXB.  However,  it  is  possible  that 
a  systematic  effect  could  be  present  during  observations  of  the 
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Figure  2.  Suzaku  HXD-PIN  Earth-occulted  data  (black  data  points)  compared 
with  the  NXB  model  spectrum  for  the  same  time  periods  of  Earth  occultation 
(red  data  points).  During  Earth  occultation,  the  only  events  should  be  due  to 
the  NXB.  Note  that  the  NXB  agrees  well  with  the  normalization  and  shape  of 
the  Earth-occulted  data  at  all  energies  >12  keV  (only  this  range  is  used  in  the 
spectral  fits  for  Coma). 

sky  that  would  not  exist  during  Earth-occulted  observations,  and 
so  it  would  not  be  included  in  the  NXB  model  or  the  estimate 
of  the  systematic  uncertainty.  Mizuno  et  al.  (2008)  attempt  to 
test  for  this  possibility  with  “blank  sky”  observations  and  find  a 
larger  effective  systematic  uncertainty.  It  is  clear  that  some  part 
of  this  uncertainty  is  due  to  the  fact  that  the  fields  are  not  entirely 
“blank”  and  that  the  source  flux  will  vary  field-to-field.  Here  we 
refer  to  relatively  bright  sources  not  considered  to  be  part  of 
the  more  uniform  CXB,  which  has  a  variance  based  on  the  XIS 
sensitivity  to  point  sources  that  can  be  taken  into  account.  When 
considering  many  observations  of  one  region  on  the  sky,  so 
that  the  contamination  from  sources  will  vary  less,  the  standard 
deviation  drops  from  5.8%  derived  from  many  fields  to  5.0%, 
both  of  which  includes  a  statistical  error  of  about  3.3%.  While  an 
additional  systematic  uncertainty,  only  in  effect  when  observing 
the  sky,  cannot  be  ruled  out,  this  drop  suggests  that  systematic 
error  estimates  derived  from  sky  observations  are  somewhat 
conservative.  Because  the  contribution  of  contaminating  sources 
to  systematic  error  estimates  is  thus  far  not  well  characterized, 
we  use  the  value  derived  from  Earth-occulted  observations 
throughout  to  avoid  overestimating  this  error.  We  rely  on  the 
assumption  that  a  full  accounting  of  contaminating  sources 
would  lead  to  a  systematic  error  estimate  similar  to  our  adopted 
value.  However,  using  the  sample  of  10  ks  exposures  of  all  blank 
sky  observations  leads  to  an  estimate  of  the  NXB  systematic 
error  of  4%,  after  subtracting  the  statistical  error  and  the 
expected  CXB  fluctuation  (see,  e.g.,  Nakazawa  et  al.  2009). 
We  consider  the  effect  of  raising  the  systematic  error  to  this 
higher  value  (for  E  <  40  keV)  in  Section  4.3;  our  results  and 
conclusions  remain  qualitatively  unchanged. 

3.2.  XMM  EPIC-pn  Spectrum 

To  constrain  the  thermal  contribution  to  the  PIN  hard 
X-ray  spectrum,  it  is  very  helpful  to  have  a  spectrum  for  the 
same  region  covered  by  the  HXD-PIN  FOV  but  extending  to 
lower  energies  where  the  thermal  emission  is  completely  domi¬ 
nant.  This  low  energy  spectrum  acts  as  a  lever  arm  on  the  thermal 
continuum  so  that  the  properties  of  the  thermal  gas  can  be  ex¬ 
tracted  simultaneously  with  a  potential  nonthermal  component. 
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Figure  3.  XMM  EPIC-pn  2-7.5  keV  wavelet  smoothed  X-ray  surface  brightness 
image  (Schuecker  et  al.  2004)  with  contours  of  constant  Suzaku  HXD-PIN 
effective  area  overlaid.  The  contours  are  spaced  at  10%  intervals  running  from 
0%  to  90%  of  the  effective  area  for  a  point  source  located  at  the  instrument 
center.  In  between  the  contours,  the  regions  labeled  Rl,  R2,  etc.,  correspond  to 
those  in  Table  1  and  in  the  text. 

Because  the  ICM  of  Coma  is  not  isothermal  and  its  projected 
temperature  varies  across  the  cluster,  a  complimentary  spectrum 
at  softer  energies  must  follow  the  spatial  sensitivity  of  Suzaku' s 
HXD.  Otherwise,  localized  regions  of  even  slightly  hotter  gas 
could  mimic  the  emission  of  a  nonthermal  source  at  hard  en¬ 
ergies.  Since  the  HXD  is  made  up  of  64  individual  collimators 
with  optical  axes  generally  aligned  to  within  4'  of  each  other, 
we  approximate  the  PIN  spatial  response  as  a  single  perfect 
collimator  with  a  total  square  FOV  of  D  =  65^5  on  a  side, 

Rc on  =  (D/2  -  6x)(D/2  -  6y)/(D/ 2)2,  (3) 

where  Rco\\  is  the  fraction  of  the  flux  detected  at  angles  of 
(0X,  6y)  from  the  optical  axis  along  the  PIN  detector  axes, 
relative  to  a  point  source  located  at  the  center  of  the  HXD 
FOV  {0X  =  0y  =0).  We  have  verified  that  this  model  fits  the 
spatial  response  of  the  PIN  very  well.  The  complimentary  XMM 
spectrum  is  constructed  based  on  this  spatial  vignetting  of  the 
PIN,  which  is  reasonable  for  our  energy  range  of  interest  (less 
than  70  keV). 

In  order  to  build  a  spectrum  that  reflects  the  PIN  vignetting 
with  good  statistics,  we  extract  spectra  from  10  regions  of 
roughly  equal  effective  area,  as  shown  in  Figure  3.  The  bound¬ 
aries  of  the  regions  are  spaced  at  intervals  of  10%  of  the  PIN 
sensitivity  to  a  central  point  source.  Because  the  solid  angle  sub¬ 
tended  by  a  region  increases  with  its  distance  from  the  cluster 
center,  it  turns  out  that  the  count  rates  of  each  of  these  XMM 
spectra  are  comparable.  The  same  response  matrix  is  used  for 
all  spectra,  epn_ef20_sdY7_medium.rsp,  and  the  auxiliary  re¬ 
sponse  files  (arfs)  for  each  region  are  generated  in  the  standard 
way  (Lumb  et  al.  2003).  The  background  spectrum  is  derived 
from  the  data  sets  compiled  by  Read  &  Ponman  (2003),  to 
which  we  also  apply  consistent  flare  cleaning  criteria.  Before 
summing  these  spectra,  weighted  by  the  average  PIN  sensitivity 
within  each  region,  we  scale  the  arfs  so  they  all  agree  with  the 
central  region  (R10)  arf  at  5  keV,  while  also  scaling  the  exposure 


Table  1 

XMM  Regions  and  Spectral  Fits 


Region 

PIN  Weight 

kT 

(keV) 

Norm.a 
(10-2  cm"5) 

R10 

0.933 

8.25  ±0.21 

1.67  ±0.01 

R9 

0.844 

8.33  ±0.14 

4.25  ±  0.02 

R8 

0.746 

8.27  ±0.15 

5.17  ±0.03 

R7 

0.647 

8.07  ±0.17 

5.09  ±  0.03 

R6 

0.547 

8.07  ±  0.23 

4.52  ±  0.04 

R5 

0.448 

7.40  ±  0.34 

3.65  ±  0.05 

R4 

0.348 

7.39  ±  0.46 

2.93  ±  0.05 

R3 

0.248 

6.99  ±  0.56 

2.84  ±  0.06 

R2 

0.147 

7.65  ±  0.44 

3.30  ±  0.05 

Rl 

0.0421 

7.45  ±  0.68 

2.76  ±  0.06 

Note. 

a  Normalization  of  the  APEC  thermal  spectrum,  which  is  given  by 
{10-14/[47t(1  +  z)2D\]}  f  neriH  dV,  where  z  is  the  redshift,  Da  is  the  an¬ 
gular  diameter  distance,  ne  is  the  electron  density,  nu  is  the  ionized  hydrogen 
density,  and  V  is  the  volume  of  the  cluster. 


times  so  the  flux  remains  unchanged.  Similarly,  normalized  arfs 
are  required  to  ensure  that  the  weighted  and  summed  arf  will 
properly  represent  the  response  of  the  final  summed  spectrum. 

We  now  describe  the  procedure  for  creating  the  summed 
observed  rate  spectrum,  summed  background  spectrum,  and 
corresponding  response.  Let  Ol-  be  the  observed  rate  spectrum 

in  spectral  channel  j  for  region  /,  and  let  Blj  be  the  corresponding 
background  spectrum.  We  constructed  the  weighted  sums 

10 

Oj  =  J2wi°j’  (4) 

i=l 

10 

A  =  Xw,'%  (5) 

i= 1 

where  wl  is  the  weight  of  region  i  in  the  PIN  spectrum,  based 
on  the  average  value  of  Equation  (3)  inside  the  region  (Table  1), 
normalized  by  the  PIN  nominal  central  point  source  sensitivity. 
Let  Rljk  be  the  response  matrix  for  region  i,  defined  such  that 
given  a  number  flux  Mlk  of  photons  per  unit  area  and  time  at 
Earth  in  spectral  channel  k ,  RljkMlk  is  the  number  of  events  per 
second  observed  in  spectral  channel  j.  In  the  nomenclature  of 
X-ray  spectral  analysis,  Rljk  is  the  “rsp”  file  for  region  i.  The 

corresponding  average  response  matrix,  Rjk,  must  be  weighted 
both  by  the  PIN  effective  area  for  the  regions  (the  weights  wt) 
and  by  the  number  of  XMM  photons  incident  on  each  area. 
To  include  the  weighting  by  the  incident  flux  on  each  area,  we 
formally  assume  that  to  first  order  the  spectra  in  all  10  regions  are 
described  by  models  M[  that  have  identical  spectral  shapes  given 
by  mk ,  but  differing  normalizations  N*.  That  is,  Mlk  =  Nlnik.  In 
our  case,  we  take  as  an  APEC  model  with  T  =  8.2  keV  and 
abundance  relative  to  solar  of  0.24,  which  is  the  best-fit  model 
to  the  XMM-Newton  EPIC-pn  spectrum  for  the  entire  Suzaku 
PIN  FOV.  While  the  temperatures  from  the  outer  five  regions 
(R5-R1 ;  Table  1)  are  somewhat  inconsistent  with  our  fiducial  T, 
the  central  regions  are  weighted  more  strongly  and  so  it  is  more 
important  to  accurately  match  their  spectral  shape  than  that  of 
the  spectra  from  R5-R1.  The  average  of  the  temperatures  from 
R10-R1  is  in  good  agreement  with  previous  measurements  of 
Coma’s  global  temperature  (Arnaud  et  al.  2001;  Watanabe  et  al. 
1999).  Also,  temperatures  in  R10,  R9,  and  R8  are  just  consistent 
at  the  90%  level  with  continuum  and  Fe  line  ratio  fits  to  the  XIS 
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data  extracted  from  those  regions,  using  the  method  described 
in  T.  Sato  et  al.  (2009,  in  preparation).  The  models  Mlk  now 
differ  only  in  overall  flux,  given  by  the  APEC  normalization  N\ 
so  each  observed  spectrum  can  be  described  as 

O) -*'=£  r)M  =  N‘  E  R>k  •  (6) 

k  k 

Similarly,  we  define  Rjk  & s 

Oj-Bj  =  N'Y^Rjkmk,  (7) 

k 

where  N  is  the  APEC  normalization  of  the  summed  spectrum. 
Substituting  Equations  (4)  and  (5)  into  Equation  (7)  yields,  after 
some  algebra, 

E  RJk,nk  =  E  mk  E  w‘  77  r>j^  (8) 

k  k  i 

so  it  is  clear  that 

10  N* 

^  =  EwV^-  (9) 

i—  1 

The  value  of  the  weighted  normalization  N  is  given  by  N  = 
uj1  Nl . 

For  all  fits  of  the  XMM  spectra,  the  energy  range  considered 
differed  slightly  from  the  nominal  2-12  keV  energy  range  due 
to  calibration  issues.  At  energies  near  ~2  keV,  there  exists  a 
sharp  edge  in  the  response  due  to  gold  in  the  mirrors,  and 
between  8  and  9  keV  there  are  variable  background  lines  due 
to  copper  and  zinc.  Neither  of  these  features  can  be  sufficiently 
accounted  for  given  the  current  calibration,  and  they  tend  to 
become  important  in  regions  with  very  good  statistics  (all  of  the 
Coma  XMM  spectra)  and/or  regions  of  low  surface  brightness 
particularly  near  the  outer  edges  of  the  detector.  We  avoid  these 
problems  by  excluding  these  features  and  only  fitting  over  the 
range  2.3  <  E  <  7.5  keV  and  9.5  <  E  <  12  keV.  Also, 
the  gain  can  vary  by  up  to  30  eV,  especially  during  the  period 
when  many  of  these  observations  took  place  (Marcus  Kirsch, 
XMM  Calibration  Document9)  and  the  redshift  found  from  fits 
differs  significantly  from  the  nominal  value  of  z  =  0.023.  We 
fit  the  gain  with  a  linear  function,  assuming  the  redshift  to  be 
0.0232,  using  the  gain  command  in  XSPEC  for  all  10  spectra. 
We  adjusted  the  gain  such  that  the  new  response  energies  E'  are 
related  to  the  original  energies  by  E'  =  E/ 1 .005 19  +  0.0103 12. 
None  of  these  calibration  issues  have  any  important  effect  on  the 
characterization  of  the  continuum  features  (IC  and  hard  thermal 
bremsstrahlung)  which  are  the  subject  of  this  paper.  However, 
not  including  these  effects  would  result  in  high  values  of  x2 
for  the  fits  due  to  the  very  good  statistics  in  the  Coma  XMM 
EPIC-pn  spectra,  and  thus  make  it  more  difficult  to  determine 
the  uncertainties  in  parameters. 

Even  after  these  adjustments,  fits  to  the  weighted  and  summed 
final  spectrum  with  any  model  produce  high  chi-squared  values 
(Xred  >  1-3).  A  close  inspection  of  the  continuum  at  various 
energies  reveals  that  the  residuals  are  slightly  larger  relative  to 
the  errors  than  would  be  expected  by  x2  statistics,  indicating 
that  we  have  underestimated  the  errors.  Because  of  the  high 
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signal  to  noise  of  the  Coma  observations,  the  statistical  errors 
no  longer  completely  dominate  over  channel-to-channel  system¬ 
atic  effects,  caused  by,  e.g.,  differing/ variable  charge  transfer 
inefficiencies  across  the  detectors  and/or  between  observations, 
which  were  obtained  over  a  two  year  period.  We  add  a  3% 
systematic  error  to  the  count  rate  of  each  channel  in  the  final 
spectrum  to  account  for  these  uncertainties,  which  leads  to  more 
reasonable  values  of  xr2  d  ~  1 . 

3.3.  XMM  EPIC-pn  and  Suzaku  HXD-PIN  Cross-Calibration 

Joint  fits  between  data  from  different  instruments /missions 
require  a  careful  consideration  of  the  relative  overall  calibration 
if  the  validity  of  fits  are  to  be  believed.  Instead  of  directly 
finding  the  cross-normalization  through  other  observations  of 
a  spectrally  simple  source,  such  as  the  Crab,  we  choose  to 
use  the  XIS0  chip  of  the  Suzaku  XRT  as  an  intermediary. 
An  advantage  of  this  method  is  that  it  does  not  require  any 
assumptions  about  the  stability  of  the  absolute  calibration 
of  each  instrument  between  calibration  observations  and  our 
observations.  To  justify  our  use  of  the  XIS0  data  to  calibrate  the 
absolute  flux  level,  we  compare  the  flux  observed  by  XIS0  to 
the  ROSAT  0.5-2  keV  flux,  which  was  derived  using  Snowden’s 
ESAS  software  package  (Snowden  &  Kuntz  1998);  these  fluxes 
agree  to  within  1%. 

The  overlapping  spatial  and  spectral  coverage  of  the  XMM 
EPIC-pn  and  XIS0  instruments  allows  a  trivial  comparison  of  the 
flux  for  a  region  on  the  sky.  We  extract  an  XIS0  spectrum  from 
the  same  region  as  XMM  spectrum  R10,  and  we  generate  rmf  and 
arf  files  for  this  region  using  the  2-7.5  keV  wavelet- smoothed 
image  created  from  the  XMM  EPIC-pn  data  (Schuecker  et  al. 
2004).  Though  the  large  XIS  point-spread  function  (~2')  will 
scatter  photons  both  into  and  out  of  this  region  to  a  much  greater 
degree  than  occurs  for  XMM ,  this  effect  is  accounted  for  in  the 
arf  and  tied  to  the  XMM  data.  So  while  spatial  inhomogeneities 
will  not  impact  the  comparison,  the  shapes  of  the  spectra  will 
not  necessarily  be  identical.  The  overall  flux,  however,  is  not 
sensitive  to  small  variations  in  the  temperature,  and  so  it  provides 
a  good  quantity  to  establish  the  XMM-X IS  cross-normalization. 
For  this  region,  we  find  that  the  XMM  flux  is  15%  below  the 
XIS0  flux  from  2-7.5  keV,  and  the  XMM  flux,  extrapolated  to 
0.5-2  keV,  is  similarly  15%  below  the  ROSAT  flux,  and  so  we 
scale  the  summed  XMM  arf  by  this  factor. 

The  cross-normalization  between  the  XIS  chips  and  the 
HXD-PIN  has  been  well  studied  for  observations  of  the  Crab 
nebula  (Ishida  et  al.  2007b).  We  adopt  their  PIN/XIS0  relative 
normalization  factor  of  1.132  zb  0.014,  increasing  the  PIN  arf, 
and  thus  lowering  the  measured  flux,  by  13.2%. 

The  associated  systematic  error  for  both  cross-normalization 
corrections  is  estimated  to  be  1-2%.  However,  the  normaliza¬ 
tion  of  the  R10  spectrum  may  differ  from  that  of  the  other  XMM 
region  spectra,  and  also  the  XIS0-PIN  relative  normalization 
may  be  different  due  to  the  fact  that  Coma  is  spatially  diffuse 
while  the  Crab  nebula  is  comparable  in  size  to  the  XIS  reso¬ 
lution.  These  issues  suggest  that  the  true  cross-normalization 
systematic  uncertainty  is  probably  larger.  We  therefore  take  the 
combined  cross-normalization  systematic  error  to  be  5%,  which 
is  about  as  large  as  can  be  reasonably  allowed  by  the  simple  con¬ 
straint  that  a  model  can  be  continuously  fit  across  the  12  keV 
boundary  between  the  XMM  and  PIN  spectra.  Specifically,  we 
vary  the  cross-normalization  until  the  average  of  the  highest 
signal-to-noise  PIN  channels,  covering  12  keV  <  E  <  16  keV, 
disagrees  with  the  model  by  ~  2a -3a. 
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3.4.  Cosmic  X-ray  Background 

We  modeled  the  CXB  spectrum  shape  following  Boldt  (1987), 
specifically  using  the  analytical  form  proposed  by  Gruber  et  al. 
(1999)  based  on  the  HEAO-1  /A2+A4  data.  This  shape  is  well 
established  over  the  energy  range  3  <  E  <  60  keV  and  has 
been  confirmed  in  subsequent  measurements  (e.g.,  Revnivtsev 
et  al.  2003;  Churazov  et  al.  2007;  Ajello  et  al.  2008).  We 
adopt  a  10%  larger  normalization  of  the  spectrum,  relative  to 
the  original  HEAO-1  determination,  to  agree  with  the  more 
recent  measurements  by  INTEGRAL  (Churazov  et  al.  2007). 
This  increase  is  further  justified  by,  and  consistent  with,  the 
(8  zb  3)%  higher  normalization  found  with  Swift  (Ajello  et  al. 
2008).  Though  these  most  recent  measurements  lie  slightly,  but 
systematically,  above  the  canonical  spectrum,  as  noted  by  Ajello 
et  al.  (2008)  they  are  not  inconsistent  with  other  observations  at 
E  >  10  keV.  At  the  peak  of  the  CXB  spectrum,  the  measurement 
precision  of  HEAO-1  is  10%  (Marshall  et  al.  1980),  and  the 
measurement  made  with  the  BeppoSAX  PDS  is  consistent  at 
the  90%  level  with  a  normalization  12%  larger  (Frontera  et  al. 
2007).  In  XSPEC,  we  model  the  CXB  as 

/  E  \“L29 

CXB(£)  =  1.056  x  1(T2  (  -j-j— ^  j 

x  e~EKAhn  keV)photons  cm'2  s"1  keV\  (10) 

where  the  normalization  is  set  by  a  2°  x  2°  solid  angle  of  the 
sky  to  match  the  provided  response  file  (see  Section  3.1). 

Cosmic  variance  due  to  large-scale  structure  depends  on  the 
solid  angle  of  the  observation  (£2  =  0.32  deg2  for  the  PIN)  and 
on  the  cut-off  flux  of  removed  point  sources  (Scut),  determined 
by  the  XMM  source  completeness  (Finoguenov  et  al.  2004)  to 
be  5^(12-70  keV)  =  2.2  x  10-13  erg  s-1  cm-2.  Since  the 
variance  <7Cxb/7cxb  oc  Q-0,5^25,  we  can  estimate  the  variance 
in  our  observation  relative  to  another  measurement  assuming  a 
log  A-log  S  relation  of  N(S )  a  S'-1,5.  Using  the  HEAO-1  A2 
estimate  (Shafer  1983;  Barcons  et  al.  2000;  Revnivtsev  et  al. 
2003)  with  Q  =  15.8  deg2,  Scat  =  8  x  10-11  erg  s-1  cm-2, 
and  ctcxb/Icxb  =  2.8%  (la),  we  find  a  variance  of  7.4%  (90% 
confidence),  which  we  take  as  an  additional  systematic  error  in 
the  PIN  flux.  To  account  for  the  10%  discrepancy  between  the 
HEAO-1  and  the  INTEGRAL  (Churazov  et  al.  2007)  and  Swift 
(Ajello  et  al.  2008)  observations,  we  also  estimate  the  standard 
deviation  of  these  measurements  to  be  7%  (90%  confidence). 
Adding  these  uncertainties  in  quadrature,  we  take  the  total 
systematic  error  in  the  CXB  normalization  to  be  10%.  Below 
20  keV,  the  CXB  emission  is  <10%  of  the  total  flux,  and  it  just 
becomes  the  dominant  source  of  emission  at  ~50  keV. 

For  the  XMM  data,  the  background  spectra  include  unresolved 
point  sources  that  make  up  the  CXB,  so  they  do  not  need  to  be 
modeled. 

3.5.  Point  Sources 

Point  sources  in  the  XMM-Newton  Coma  mosaic  have  already 
been  identified  by  Finoguenov  et  al.  (2004),  who  also  give 
their  count  rates  in  three  energy  bands  (0.5-1  keV,  1-2  keV, 
and  2-4.5  keV).  For  each  of  the  72  sources,  we  assume  the 
spectrum  to  be  described  by  an  unabsorbed  power  law  and  fit 
this  model  to  each  spectrum.  We  found  that  the  sum  of  all  these 
models,  weighted  by  wl  according  to  their  positions,  could  be 
more  concisely  described  by  the  sum  of  two  power  laws  with 
photon  indices  2.1  and  1.6  and  normalizations  8.54  x  10  5  and 
1.23  x  10-4  photons  cm-2  s-1  keV-1  at  1  keV,  respectively. 


While  a  simple  power-law  description  poorly  characterizes 
some  of  the  sources,  care  is  taken  to  ensure  that  individual 
fits,  when  extrapolated  to  high  energies,  are  not  unphysical. 
Their  composite  spectrum  accounts  for  <1%  of  the  XMM  flux 
and  is  therefore  unimportant  relative  to  other  systematic  effects. 
For  this  reason  we  do  not  go  to  the  extra  effort  to  exclude  the 
sources  from  the  XMM  spectra.  Assuming  the  spectral  fits  are 
reasonably  valid,  the  point  sources  account  for  ~10%  of  the 
CXB  at  PIN  energies.  During  fits  of  the  joint  XMM  and  Suzaku 
data,  we  include  this  point  source  composite  model  for  both 
spectra. 

The  brightest  of  these  point  sources  is  X  Comae,  a  background 
AGN  with  a  flux  ~10  times  brighter  than  any  other  source  in 
the  field.  From  XMM  RGS  observations,  it  is  known  to  have  a 
steep  spectrum  (r  ~  2.4)  and  to  vary  in  flux  by  about  a  factor 
of  2  over  the  course  of  one  year  (Takei  et  al.  2007).  However, 
due  to  its  position,  nearly  90%  of  the  flux  from  X  Comae  is 
not  detected  by  the  HXD,  so  this  source  does  not  significantly 
contribute  to  the  PIN  spectrum. 

4.  SPECTRAL  FITS 

In  our  spectral  fits,  all  model  components  are  absorbed  by 
the  neutral  hydrogen  column  density  toward  Coma,  NH  = 
9  x  1019  cm2  (average  of  values  derived  from  Dickey  & 
Lockman  1990  and  Kalberla  et  al.  2005),  though  this  absorption 
is  negligible  at  energies  above  2  keV.  In  general,  we  characterize 
the  dominant  thermal  emission  in  XSPEC  with  the  APEC  model 
for  E  <  40  keV  and  with  the  MeKa  model  for  E  >  40  keV.  In 
the  correct  version  of  XSPEC,  the  APEC  and  MeKaL  models  are 
undefined  above  50  keV,  though  the  MeKa  and  bremsstrahlung 
models  are  defined.  We  tie  the  parameters  of  the  MeKa  model 
to  the  APEC  parameters,  except  for  the  MeKa  normalization, 
which  we  reduce  relative  to  the  APEC  normalization  by  5.5% 
to  bring  the  models  into  agreement  at  high  energies.  Also, 
the  abundances  relative  to  solar  and  the  redshift  are  fixed, 
to  0.24  and  0.0232  respectively  (see  Section  3.2).  This  value 
for  the  abundance  is  based  on  fits  to  the  final  weighted  and 
summed  XMM  spectrum  alone,  and  the  best-fit  abundances  of 
all  the  individual  spectra  from  the  10  regions  is  also  consistent 
with  this  value.  Since  we  are  interested  in  continuum  features, 
the  exact  choice  for  the  abundance  does  not  strongly  affect 
the  results.  The  spectral  fitting  results  are  summarized  in 
Table  2. 

4.1.  Joint  XMM-Newton  and  Suzaku  Spectral  Fits  Without 
Considering  Systematic  Errors 

We  simultaneously  fit  the  Suzaku  HXD-PIN  and  XMM- 
Newton  EPIC-pn  spectra  for  the  PIN  FOV.  First,  we  consider 
only  a  single-temperature  (single- T)  fit,  in  order  to  establish 
whether  the  addition  of  a  nonthermal  component  actually 
improves  the  fit  (Figure  4).  We  find  a  best-fit  temperature  of 
8.45  ±  0.06  keV,  which  is  in  general  agreement  with  similar  fits 
to  the  PIN  data  (8 . 3  =b  0. 3  keV)  and  XMM  data  (8 . 37  ±  0. 1 2  keV) 
individually.  Note  that  the  dip  at  15  keV  is  a  known  problem 
with  the  current  NXB  model  (Mizuno  et  al.  2008).  Since 
each  spectrum  is  individually  described  by  the  same  average 
temperature,  the  existence  of  excess  emission  at  hard  energies 
is  unlikely.  While  all  of  these  temperatures  are  slightly  higher 
than  the  cluster- wide  average  temperature  of  8.2  keV  (Hughes 
et  al.  1993),  the  energy  range  in  this  and  similar  fits  typically 
extends  to  energies  below  2  keV  and  thus  includes  more  low 
temperature  gas. 
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Table  2 

Fits  to  Joint  XMM- PIN  Spectra 


Model 

kT 

(keV) 

Norm.a 

(cm-5) 

Tor  kT 
(b) 

Norm.c 

X2/dof 

Single  T 

8.45  ±  0.06 

0.218  ±0.001 

1676.05/1689 

T+ ICd 

8.42  ±  0.06 

0.218  ±0.001 

-1.6 

(4.6  ±  3.5)  x  10“9 

1671.29/1688 

T+  ICd 

8.45  ±  0.07 

0.217  ±  0.002 

2.0 

(2+22)  x  10“4 

1676.18/1688 

IT 

8.0 

0.17 

10.1 

0.05 

1672.34/1687 

^map 

1684.35/1690 

Notes. 

a  See  the  note  following  Table  1 . 

b  Value  is  T  for  the  T  +IC  model  and  kT  (in  keV)  for  the  2 -T  model. 

c  Value  is  the  normalization  of  the  power-law  component  for  the  T  +  IC  model,  which  is  the  photon  flux  at  a  photon  energy  of  1  keV  in  units  of  photons 
cm-2  s-1  keV-1.  For  the  2 -T  model,  the  value  is  the  normalization  of  the  second  APEC  thermal  model  (see  the  note  following  Table  1)  in  units  of 
cm-5. 

d  Value  of  T  is  fixed  when  deriving  errors. 
e  Parameters  unconstrained. 
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Figure  4.  Suzaku  HXD-PIN  spectrum  (E  >  12  keV)  and  the  combined  XMM 
spectrum  (E  <  12  keV)  corresponding  to  the  spatial  sensitivity  of  the  PIN. 
Shown  as  solid  lines  are  the  best-fit  models  for  a  single-  T  thermal  component. 
The  thermal  model  (“APEC,”  green)  is  nearly  coincident  with  the  data,  though 
falling  below  it  at  higher  energies.  Also  included  for  all  joint  fits  are  the  the  total 
spectrum  for  the  “XMM  Point  Sources”  (red)  and  the  Cosmic  X-ray  Background 
(“CXB,”  purple),  the  latter  of  which  only  applies  to  the  PIN  spectrum  since  the 
CXB  is  subtracted  from  the  XMM  data  along  with  the  NXB. 

(A  color  version  of  this  figure  is  available  in  the  online  journal.) 


The  addition  of  a  power-law  nonthermal  component  produces 
a  formally  better  description  of  the  spectra,  according  to  the  f- 
test,  improving  the  overall  fit  (Table  2),  but  only  for  a  photon 
index  T  <  0.  Allowing  the  temperature  and  power-law  photon 
index  to  vary  along  with  each  component’s  normalization,  we 
find  T  =  8.42  zb  0.06  keV  and  T  =  —1.6,  though  T  is 
poorly  constrained.  If  we  fix  T  to  this  best-fit  value,  the  IC 
component  is  significant  at  the  2.2cr  level  without  considering 
the  effect  of  systematic  uncertainties.  However,  this  photon 
index  is  completely  inconsistent  with  the  spectral  index  of  the 
radio  halo  (T  >  1.8,  Giovannini  et  al.  1993).  While  we  might 
expect  a  flatter  spectrum  for  IC  emission,  since  the  hard  X-ray 
photons  are  emitted  by  somewhat  lower  energy  electrons  than 
the  radio  emission,  and  the  radio  spectrum  flattens  at  lower 
frequencies  (Thierbach  et  al.  2003),  a  rising  IC  spectrum  with 
energy  is  completely  unexpected  and  unphysical.  The  power- 
law  fit,  in  contrast  to  finding  an  actual  power-law  signature  in 
the  data,  is  instead  compensating  for  a  slight  excess  at  high 


Energy  (keV) 

Figure  5.  Suzaku  HXD-PIN  spectrum  (E  >  12  keV)  and  the  combined  XMM 
spectrum  (E  <  12  keV)  corresponding  to  the  spatial  sensitivity  of  the  PIN. 
Shown  as  solid  lines  are  the  best-fit  models  for  a  single- T  thermal  component 
plus  a  nonthermal  component.  The  thermal  model  (“APEC,”  green)  is  nearly 
coincident  with  the  data,  though  falling  below  it  at  higher  energies.  The 
nonthermal  model  (“Power  Law,”  light  blue)  is  the  faintest  model  component 
for  both  spectra,  and  the  photon  index  is  fixed  at  T  =  2.0.  The  other  two 
components  are  described  in  Figure  4. 

(A  color  version  of  this  figure  is  available  in  the  online  journal.) 

energies  while  minimizing  its  impact  on  the  overall  fit  at  lower 
energies.  Notice  that  the  residuals  in  Figure  4  above  40  keV 
lie  systematically,  if  not  significantly,  above  the  model.  This 
excess  at  energies  above  40  keV  can  be  explained  as  a  ~2% 
underestimate  of  the  NXB,  as  suggested  by  the  Earth-occulted 
spectrum  (see  Section  3.1).  Increasing  the  background  level  by 
2%  for  E  >  40  keV  results  in  a  best-fit  power-law  component 
very  similar  to  the  model  used  for  the  XMM  point  sources,  with 
r  ~  1.6,  but  it  is  not  significant  at  even  the  la  level.  A  similar 
result  is  found  if  the  photon  index  is  fixed  at  F  =  2  and  the  NXB 
above  40  keV  is  not  increased;  this  fit  is  shown  in  Figure  5.  In 
this  case,  the  fit  is  not  improved  by  the  addition  of  a  power-law 
component  to  the  model. 

Interestingly,  a  two-temperature  (2-7)  model  for  the  ICM 
yields  only  a  slightly  better  fit  to  the  data  than  does  the 
single-T  model  (see  Figure  6),  though  the  addition  of  a  second 
temperature  component  is  probably  not  formally  justified.  This 
result  is  mainly  due  to  the  addition  of  the  3%  systematic  error  to 
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Figure  6.  Suzaku  HXD-PIN  spectrum  (E  >  12  keV)  and  the  combined  XMM 
spectrum  (E  <  12  keV)  corresponding  to  the  spatial  sensitivity  of  the  PIN. 
Shown  as  solid  lines  are  the  best-fit  models  for  a  2-T  thermal  component. 
The  thermal  model  (“APEC+APEC,”  green)  is  nearly  coincident  with  the  data, 
though  falling  below  it  at  higher  energies.  The  other  two  components  are 
described  in  Figure  4. 


the  XMM  spectrum.  Without  including  that  error,  a  2-T  model 
produces  a  clearly  improved  fit  over  a  single-T model,  indicating 
that  the  addition  of  this  error  is  somewhat  obscuring  evidence 
for  a  multitemperature  continuum.  In  either  case,  the  two 
temperatures  are  not  strongly  constrained,  but  they  are  broadly 
consistent  with  the  spatial  variations  in  Coma’s  temperature  (see 
Section  4.2  and  Figure  7).  Therefore,  even  before  systematic 
errors  are  considered,  the  case  for  the  inclusion  of  a  nonthermal 
component  is  not  strongly  motivated. 


4.2.  Multiple  Thermal  Components  in  Coma 

In  most  previous  attempts  to  measure  a  nonthermal  compo¬ 
nent  in  the  hard  spectrum  of  Coma,  the  thermal  emission  was 
modeled  as  a  single- T plasma  characteristic  of  the  average  global 
state  of  the  ICM.  However,  Coma  is  known  to  host  temperature 
variations  (e.g.,  Honda  et  al.  1998).  Generally,  all  clusters  ex¬ 
hibit  a  multitemperature  ICM  (Cavagnolo  et  al.  2008),  and  this 
is  especially  true  of  merging  clusters  like  Coma,  which  tend  to 
host  hot  regions  due  to  shocks  (Markevitch  et  al.  1998).  At  hard 
energies,  where  the  exponential  turnover  in  the  bremsstrahlung 
continuum  is  especially  well  sampled,  even  weak  higher  tem¬ 
perature  components  can  significantly  contribute  to  the  flux. 
Also,  these  components  would  lead  to  a  higher  average  temper¬ 
ature  for  the  ICM  than  if  the  average  cluster  temperature  were 
determined  from  the  spectrum  at  softer  energies,  such  as  from 
0.5-10  keV. 

In  the  previous  section,  we  found  that  a  two  tempera¬ 
ture  model  provided  a  slightly  better  description  than  did  a 
single- T  model  of  joint  fits  to  the  XMM-Newton  and  Suzaku 
data,  especially  when  ignoring  the  3%  systematic  error  applied 
to  the  XMM  data.  This  may  indicate  that  there  are  multiple 
temperature  components  in  Coma.  The  multiple  components 
could  occur  along  the  line  of  sight,  or  in  the  plane  of  the 
sky,  or  locally  (the  gas  might  be  multiphase).  In  fact,  previous 
temperature  maps  show  that  Coma  certainly  has  temperature 
structure  which  is  likely  associated  with  mergers  (Watanabe 
et  al.  1999).  Here,  we  test  whether  this  temperature  struc¬ 
ture  alone  could  reproduce  the  observed  Suzaku  PIN  spectrum 
of  Coma,  without  any  nonthermal  emission.  From  the  XMM- 
Newton  EPIC-pn  mosaic  of  Coma,  we  constructed  a  tempera¬ 
ture  map  on  a  16  x  16  grid  with  cell  size  of  4'3  on  a  side.  Each 
of  the  spectra  were  fit  with  a  single- T  APEC  model  to  produce 
a  temperature  map  that  covers  the  Coma  mosaic,  as  shown  in 
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Figure  7.  XMM-Newton  temperature  map  across  Coma  with  HXD-PIN  contours  of  constant  PIN  effective  area  overlaid  at  10%  intervals.  The  XMM-Newton  spectra 
were  fit  in  square  spatial  regions  4f3  on  a  side.  The  temperatures,  given  in  keV  by  the  color  bar,  are  accurate  to  either  a  few  tenths  of  a  keV  (in  the  center)  or  1-2  keV 
in  lower  surface  brightness  regions.  Temperatures  shown  here  were  determined  from  fits  to  the  0.5-14  keV  spectrum  in  each  region. 
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Figure  8.  Suzaku  HXD-PIN  spectrum  (E  >  12  keV)  and  the  combined  XMM 
spectrum  (E  <  12  keV)  corresponding  to  the  spatial  sensitivity  of  the  PIN. 
Shown  as  solid  lines  are  the  combined  spectra  of  the  best-fit  models  from  the 
temperature  map.  The  thermal  model  (“Tmap,”  green)  is  nearly  coincident  with 
the  data,  though  falling  below  it  at  higher  energies.  The  other  two  components 
are  described  in  Figure  4. 

(A  color  version  of  this  figure  is  available  in  the  online  journal.) 

Figure  7.  We  weighted  these  model  fits  by  the  PIN  spatial  sensi¬ 
tivity  and  combined  them.  This  resulting  model  was  compared 
with  the  PIN  spectrum  (Table  2,  row  labeled  4Tmap”).  Note  that 
only  the  overall  normalization  of  the  rmap  model  was  allowed 
to  be  fit,  to  compensate  for  a  loss  of  flux  due  to  incomplete 
coverage  of  the  map  across  the  HXD  FOV.  Also,  the  spectral 
shape  and  normalization  of  each  of  the  thermal  models  was 
the  same  as  given  by  the  XMM-Newton  temperature,  and  each 
model  was  simply  weighted  by  the  average  PIN  effective  area 
at  that  position. 

This  rmap  model  provides  a  good  fit  to  the  PIN  spectrum 
with  no  adjustable  parameters  that  affect  the  spectral  shape 
(Figure  8).  While  the  value  of  x  2  is  worse  than  for  all  previous 
fits,  it  is  only  slightly  worse,  which  is  not  unexpected  since  it  is 
not  really  a  fit  at  all.  The  fact  that  XMM-derived  temperatures 
extrapolated  to  PIN  energies  are  sufficient  to  fully  account  for 
the  PIN  spectrum  further  suggests  that  no  nonthermal  hard  X- 
ray  excess  has  been  detected  with  the  PIN,  especially  below 
40  keV  and  probably  below  70  keV.  Also,  simulating  spectra 
of  similar  quality  to  our  data,  assuming  the  rmap  model  for  the 
underlying  source,  yields  a  joint  single-T  fit  consistent  with  that 
found  from  the  actual  data,  with  T  =  8.51  ±  0.06  keV. 

4.3.  Systematic  Errors  in  the  Spectrum 

We  explicitly  consider  the  systematic  error  for  three  quanti¬ 
ties:  the  PIN  non-X-ray  background,  the  XMM-Suzaku  cross¬ 
normalization  factor,  and  the  normalization  of  the  CXB  as  mod¬ 
eled  for  the  PIN  spectrum.  To  test  the  effect  of  these  system¬ 
atic  errors,  we  vary  the  relevant  quantity  up  and  down  by  our 
estimate  of  the  90%  systematic  error,  and  evaluate  the  resulting 
change  in  best-fit  model  parameters.  The  detection  of  a  nonther¬ 
mal  component  cannot  be  claimed  unless  it  remains  robust  to 
variations  of  these  quantities  within  their  systematic  errors.  Be¬ 
cause  the  largest  error  is  in  the  normalization  of  the  PIN  NXB, 
we  first  increase  it  by  2.3%  from  12-40  keV  and  4%  from  40- 
70  keV  and  repeat  the  single- T  plus  nonthermal  model  fit.  The 
new  best-fit  IC  normalization  is  pushed  to  zero.  Even  before 
considering  the  effect  of  other  systematic  errors,  from  this  exer¬ 
cise  alone  it  is  clear  that  we  do  not  detect  nonthermal  emission 
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Figure  9.  Suzaku  HXD-PIN  spectrum  {E  >  12  keV)  and  the  combined  XMM 
spectrum  (E  <  12  keV)  corresponding  to  the  spatial  sensitivity  of  the  PIN. 
The  PIN  NXB  is  increased  to  its  90%  confidence  limit,  which  is  2.3%  for 
E  <  40  keV  and  4%  for  E  >  40  keV.  The  thermal  model  (“APEC,”  green)  is 
nearly  coincident  with  the  data,  though  falling  below  it  at  higher  energies.  The 
other  two  components  are  described  in  Figure  4.  Note  that  the  residuals  above 
the  model  that  exist  in  the  previous  fits  at  E  >  40  keV  have  disappeared. 

(A  color  version  of  this  figure  is  available  in  the  online  journal.) 

in  the  HXD-PIN  spectrum,  given  the  current  uncertainty  in  the 
NXB  normalization.  This  fit,  with  a  temperature  for  the  thermal 
component  of  8.33  =b  0.06  keV,  is  shown  in  Figure  9.  In  fact,  the 
value  of  x2  is  slightly  lower  (x2  =  1672.25  for  1689  dof)  than 
for  the  nominal  PIN  NXB  single- T  fit.  Notice  that  the  residu¬ 
als  above  the  model  for  E  >  40  keV,  seen  in  all  the  previous 
spectral  fits,  have  disappeared. 

Though  we  cannot  claim  to  detect  nonthermal  emission,  we 
can  derive  an  upper  limit  to  its  flux  based  on  joint  fits  to 
the  spectra,  including  systematic  errors  in  the  following  way. 
First,  for  an  assumed  photon  index  which  we  fix,  we  find 
the  nominal  normalization  Anom  and  corresponding  90%  upper 
bound  A^nom  °f  ^e  nonthermal  component  for  a  single- T  plus 
power-law  model,  allowing  the  temperature  and  normalization 
of  the  thermal  component  to  vary.  Then,  for  each  systematically 
uncertain  quantity,  we  set  that  quantity  to  the  limit  bounded  by 
the  systematic  error  in  the  sense  that  increases  the  value  of  the 
nonthermal  normalization  AsyM-,  and  we  fit  for  it  and  its  90% 
upper  bound  N^s  • .  The  statistical  and  systematic  errors  of  the 
power-law  normalization  are  then  given  by 

8Stat  =  KL-Naom,  (11) 

and 

8sys,i  =  <s,i  -  Kom  =  Nsys,i  -  Vnom,  (12) 

respectively.  The  final  90%  upper  limit  is  then  given  by 


<  =  Vnom  +  at-  (13> 

We  add  each  systematic  error  contribution  in  quadrature  because 
it  is  unlikely  that  we  chose  normalizations  for  these  three  quan¬ 
tities  such  that  each  one  disfavors  the  detection  of  nonthermal 
emission  in  the  most  severely  possible  way. 

The  upper  limits  for  a  range  of  assumed  photon  indices  is 
provided  in  Table  3,  and  in  Figure  10  we  show,  for  T  =  2,  the 
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Table  3 

90%  Upper  Limits  on  IC  Flux 


r 

Norm.a 

(10-3  photons  keV-1  cm-2  s-1) 

Flux  (12-70  keV) 

(10-12  erg  s-1  cm-2) 

Flux  (20-80  keV) 

(10-12  erg  s-1  cm-2) 

BeppoSAX  Detection13 
(10-12  erg  s-1  cm-2) 

1.0 

0.155 

14.4 

14.9 

l.i 

0.220 

14.2 

14.4 

1.2 

0.311 

14.1 

13.9 

1.3 

0.439 

13.9 

13.4 

1.4 

0.617 

13.8 

12.9 

1.5 

0.860 

13.5 

12.3 

1.6 

1.18 

13.1 

11.6 

1.7 

1.58 

12.4 

10.7 

1.8 

2.04 

11.4 

9.51 
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Notes. 

a  Normalization  of  the  power-law  component  for  the  T  +  IC  model,  which  is  the  photon  flux  at  a  photon  energy  of  1  keV. 
b  Flux  (20-80  keV),  as  reported  in  Fusco-Femiano  et  al.  (2004). 
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Figure  10.  Suzaku  HXD-PIN  spectrum  (L  >12  keV)  and  the  combined  XMM 
spectrum  (E  <  12  keV)  corresponding  to  the  spatial  sensitivity  of  the  PIN.  All 
quantities  with  systematic  uncertainties  (PIN  CXB  and  NXB,  the  XMM-Suzaku 
cross  calibration)  are  set  to  their  90%  confidence  limit  in  the  direction  that  favors 
the  addition  of  a  power-law  model  component  to  describe  the  data.  The  thermal 
model  (“APEC,”  green)  is  nearly  coincident  with  the  data,  though  falling  below 
it  at  higher  energies.  The  nonthermal  model  (“Power  Law,”  light  blue)  is  shown 
for  T  =  2  at  its  90%  confidence  upper  limit  value.  The  other  two  components 
are  described  in  Figure  4. 

(A  color  version  of  this  figure  is  available  in  the  online  journal.) 

resulting  best  fit  with  all  three  systematic  quantities  set  at  the 
limit  of  their  90%  confidence  range.  Over  the  PIN  energy  band 
(12-70  keV),  the  flux  is  relatively  independent  of  photon  index. 
To  compare  our  results  to  the  most  recent  previous  detections  of 
nonthermal  emission  in  Coma,  we  also  give  the  upper  limit  on 
the  nonthermal  flux  in  the  20-80  keV  band,  which  is  6.0  x  10-12 
erg  s-1  cm-2  for  T  =  2.  This  limit  is  inconsistent  with  the  RXTE 
(Rephaeli  &  Gruber  2002)  and  BeppoSAX  (Fusco-Femiano  et  al. 
2004)  detections  by  about  a  factor  of  2.5,  but  at  the  same  level 
as  the  upper  limit  derived  by  Rossetti  &  Molendi  (2004)  from 
the  BeppoSAX  data.  If  we  adopt  a  4%  systematic  error  for  the 
PIN  NXB  instead  of  2.3%  for  E  <40  keV,  which  would  better 
agree  with  that  derived  from  “blank  sky”  observations,  then  our 


upper  limit  increases  by  35%.  Similarly,  if  we  also  increase 
the  CXB  and  XMM-Suzaku  cross-normalization  to  18%  and 
10%,  respectively,  our  upper  limit  for  T  =  2  would  increase 
by  50%.  In  either  case,  our  upper  limit  still  excludes  the  RXTE 
and  BeppoSAX  detections,  if  FOV  differences  are  ignored  (see 
Section  5  for  a  more  meaningful  comparison). 

5.  IMPLICATIONS  AND  DISCUSSION 

After  modeling  all  the  known  possible  contributions  to  the 
2-70  keV  spectrum,  simultaneously  fitting  for  the  parameters 
of  thermal  and  nonthermal  spectral  components,  and  taking 
into  account  the  systematic  uncertainty  of  the  PIN  NXB,  we 
do  not  see  evidence  for  IC  emission  in  Coma  at  our  level  of 
sensitivity.  We  therefore  derive  an  upper  limit  to  nonthermal, 
hard  X-ray  emission  through  a  careful  consideration  of  the 
maximum  effect  of  systematic  uncertainties  on  our  ability  to 
detect  a  nonthermal  signal.  This  conservative  upper  limit  is 
similar  to  that  derived  by  Rossetti  &  Molendi  (2004)  and  is 
inconsistent  with  claimed  detections  using  RXTE  (Rephaeli  & 
Gruber  2002)  and  BeppoSAX  (Fusco-Femiano  et  al.  2004)  by 
approximately  a  factor  of  2.5.  However,  it  should  be  noted  that 
we  do  not  include  potentially  lost  emission  due  to  PIN  vignetting 
from  any  of  our  flux  upper  limits  relative  to  the  larger  FOVs  of 
RXTE  and  BeppoSAX ,  which  have  collimator  FWHM  of  1°  and 
1?3,  respectively. 

If  IC  emission  follows  the  radio  synchrotron  emission  (as 
derived  from  the  point  source- subtracted  radio  image  from  Deiss 
et  al.  1997),  as  it  would  for  a  uniform  B  field  throughout 
the  cluster,  our  upper  limits  imply  a  total  IC  flux  1.7-2x 
larger  would  be  found  inside  an  RXTEIBeppoSAX-likQ  FOV. 
We  also  consider  a  more  detailed  spatial  distribution  for  the 
underlying  IC  emission,  derived  from  the  reacceleration  model 
of  Brunetti  &  Blasi  (2005),  in  which  relativistic  protons  collide 
with  electrons  in  the  ICM  that  are  then  reaccelerated  by  Alfven 
waves  due  to  cluster  mergers.  Given  the  radio  spectrum  of  Coma, 
this  model  predicts  that  the  smaller  FOV  of  the  Suzaku  HXD- 
PIN  would  lead  to  an  underestimate  of  the  nonthermal  flux  by  a 
factor  of  2-2.5  (possibly  three  under  extreme  circumstances).  If 
this  model  for  the  spatial  distribution  of  the  nonthermal  emission 
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is  correct,  then  our  upper  limit  is  just  consistent  with  the 
measurements  of  Rephaeli  &  Gruber  (2002)  and  Fusco-Femiano 
et  al.  (2004;  however,  the  Brunetti  &  Blasi  2005  model  actually 
predicts  a  nonthermal  flux  considerably  below  the  BeppoSAX 
and  RXTE  detections).  Because  any  spatial  variation  of  the 
magnetic  field  strength  is  unknown,  a  direct  comparison  between 
these  missions  is  not  possible  with  any  precision.  Under  the 
reasonable  assumption  that  B  decreases  with  radius,  our  upper 
limit  will  be  >2  times  larger,  so  our  result  cannot  definitively 
rule  out  the  detections  discussed  above.  Regardless  of  this 
issue,  the  upper  limit  is  approximately  the  same  as  or  slightly 
higher  than  the  upper  limit  range  found  by  Rossetti  &  Molendi 
(2004). 

However,  the  BeppoSAX  PDS  measures  a  20-80  keV  flux 
for  the  Crab  of  1.23  x  10-8  erg  s-1  cm-2  (Kirsch  et  al.  2005), 
while  the  Suzaku  PIN  flux  over  this  energy  range  is  1.56  x 
10-8  erg  s-1  cm-2,  after  applying  the  13.2%  correction  to  bring 
the  PIN  spectrum  into  agreement  with  the  XIS  fluxes  (Ishida 
et  al.  2007b).  This  21%  flux  difference  implies  our  upper  limit 
would  be  4.7  x  10-12  erg  s-1  cm-2  on  the  BeppoSAX  scale, 
which  is  on  the  lower  end  of  the  range  estimated  by  Rossetti  & 
Molendi  (2004).  Also,  even  if  we  only  detect  one-third  of  the 
total  emission  observed  by  the  BeppoSAX  PDS,  we  would  just 
barely  exclude  the  nominal  value  of  the  Fusco-Femiano  et  al. 
(2004)  measurement. 

Assuming  the  differing  measurements  of  nonthermal  emis¬ 
sion  are  not  due  to  the  IC  radiation  having  a  larger  extent, 
what  might  be  the  cause  of  this  discrepency?  While  it  could  be 
explained  by  a  greater  point  source  contamination  at  hard  ener¬ 
gies  for  the  RXTE  and  BeppoSAX  missions  due  to  their  larger 
FOVs,  most  likely  we  differ  in  our  results  because  of  different 
considerations  of  the  thermal  gas.  For  both  detections,  the  gas 
temperature  was  found  to  be  lower  than  our  nominal  value  of 
8.45  keV.  Fixing  the  gas  temperature  to  their  assumed  values 
in  our  fits  yields  a  T  =  2.0  nonthermal  component  signifi¬ 
cance  greater  than  4a,  without  including  systematic  effects,  for 
T  =  7.67  keV  (RXTE)  and  T  =  8.2  keV  Q BeppoSAX );  how¬ 
ever,  these  fits  are  poor  relative  to  fits  in  which  the  temperature 
is  a  free  parameter.  Though  these  temperatures  differ  from  our 
best-fit  value  by  only  a  few  percent,  the  exponential  decline 
of  bremsstrahlung  continuum  at  high  energies  amplifies  even 
small  differences.  The  lower  measurements  of  the  ICM  temper¬ 
ature  appear  not  to  be  due  to  the  inclusion  of  data  at  low  energies 
(E  <  1  keV),  which  can  bias  average  temperature  estimates  low. 
Most  likely,  the  larger  FOVs  of  RXTE  and  BeppoSAX  allowed 
the  inclusion  of  emission  from  more  cool  gas  in  the  cluster  out¬ 
skirts  than  was  observed  by  Suzaku.  This  emission  would  serve 
to  lower  the  average  observed  temperature,  which  is  primarily 
determined  from  emission  at  lower  energies  (E  <  10  keV).  But, 
as  evidenced  by  the  temperature  map  in  Figure  7,  a  distribution 
of  higher-than-average  temperature  regions  can  effectively  in¬ 
crease  the  average  gas  temperature  observed  at  high  energies, 
as  first  seen  by  Neumann  et  al.  (2003). 

We  take  the  reasonably  good  agreement  between  the  thermal 
models  derived  from  the  XMM  temperature  map  and  the  PIN 
spectrum  to  mean  that  we  essentially  only  detect  thermal 
emission  from  Coma  out  to  70  keV.  This  result  is  fully  consistent 
with  recent  INTEGRAL  detections  of  extended  hard  X-ray 
emission.  Renaud  et  al.  (2006)  performed  a  point-by-point 
spectral  comparison  between  XMM -derived  and  INTEGRAL- 
derived  temperatures  and  found  that  they  followed  a  strict  one- 
to-one  correlation.  Similarly,  Eckert  et  al.  (2007)  characterized 
a  surface  brightness  excess  relative  to  the  XMM  data,  which  they 


found  to  be  best  described  by  extended  hot,  thermal  emission 
at  a  T  ~  12  =b  2  keV.  This  excess  coincides  with  the  hotter 
temperatures  (T  ~  10-11  keV)  to  the  west  of  the  PIN  pointing 
center  in  Figure  7. 

From  our  upper  limit  on  the  flux  of  IC  emission,  we  can  derive 
a  lower  limit  on  the  average  magnetic  field  strength  B  as  shown 
by  Harris  &  Romanishin  (1974).  Equation  (2)  refers  to  the  total 
energy  emitted  for  one  electron;  it  is  more  useful  to  consider  the 
ratio  of  monochromatic  fluxes  ER(vR),  Fx(vx)  at  frequencies 
vR,  vx,  for  a  power-law  distribution  of  electrons,  from  which 
we  can  derive  an  expression  for  the  magnetic  field  as 


/  JZ  \  2/(p+l)  / 

B  =  C(p)(l  +  j)(p+5V(/,+1)  |  _  |  (  VJL 

\FX  /  W 


(p-1)/(p+1) 

(14) 


where  p  is  the  index  of  the  electron  distribution  N(E)  oc  E~p 
and  is  related  to  the  spectral  index  a  (Fv  oc  v~a)  by  p  = 
2a  +  1.  The  value  of  the  proportionality  constant  C(p )  can  be 
found  from  the  ratio  of  the  synchrotron  flux  (Fongair  1994, 
Equation  (18.49))  to  the  IC  flux  (Rybicki  &  Fightman  1979, 
Equation  (7.31)).  Assuming  that  the  electron  energy  distribution 
does  not  turn  over  significantly  at  low  energies  and  that  a  =  1, 
we  find  B  >  0.15  /i G.  This  limit  is  still  below  the  equipartition 
value  of  0.5  juG  (Giovannini  et  al.  1993),  and  it  is  well  below 
the  line-of-sight  estimates  of  a  few  p G  derived  from  Faraday 
rotation  measure  (RM)  studies  (Feretti  et  al.  1995).  Note  that 
the  Faraday  RM  magnetic  field  estimates  are  sensitive  to  the 
B  field  geometry  and  may  imply  a  field  strength  larger  than 
the  volume-averaged  value  if  B  is  preferentially  aligned  along 
filamentary  structures  on  small  scales  (Petrosian  2001).  Also,  we 
are  unable  to  put  interesting  constraints  on  the  relativistic  energy 
budget  of  Coma,  since  our  lower  limit  includes  the  equipartition 
estimate  of  B,  which  defines  the  minimum  energy  in  relativistic 
components  and  would  not  imply  a  significant  contribution  to 
the  energy  budget  of  Coma. 

It  has  been  noted  that  the  hard  excess  detected  by  Eckert 
et  al.  (2007)  also  corresponds  to  the  peak  in  the  point  source- 
subtracted  image  from  Deiss  et  al.  (1997),  potentially  indicating 
that  the  hard  emission  could  in  fact  be  nonthermal  in  origin 
(Eckert  et  al.  2008b).  In  fact,  we  suspect  that  this  peak,  which 
appears  tantalizingly  like  a  small  radio  relic,  is  not  a  true 
feature  of  the  halo,  but  instead  is  the  result  of  imperfect  source 
subtraction.  Due  to  the  large  beam  size  used  to  create  the  diffuse 
radio  image,  extended  emission  from  radio  galaxies  might  not 
have  been  properly  subtracted  using  a  point  source  list.  We 
point  out  that  the  strongest  radio  source  (1256+282  or  5C  4.81, 
centered  on  NGC  4869)  in  Coma  is  near  this  position,  is  a  head- 
tail  radio  galaxy  with  a  steeper  spectral  index  than  rest  of  the 
halo  (Giovannini  et  al.  1993),  and  that  its  tail  extends  to  the 
west  (O’Dea  &  Owen  1985,  see  Figure  l(j))  and  turns  north 
(Venturi  et  al.  1990,  see  Figure  2).  Subtracting  a  point  source 
from  this  morphology  would  leave  a  residual  very  similar  to 
that  in  the  Deiss  image.  Therefore,  any  relation  between  the 
location  of  hard  emission  and  this  radio  feature  probably  should 
be  regarded  as  coincidence. 

As  the  calibration  of  the  NXB  model  improves,  constraints 
on  IC  emission  in  Coma  will  tighten,  possibly  leading  to  a 
detection.  The  uncertainty  in  the  current  background  model, 
“bgd_d,”  is  more  than  a  factor  of  2  lower  than  the  original 
model.  However,  the  existence  of  nonthermal  emission  in  the 
Coma  cluster  may  have  to  be  determined  by  future  missions — 
particularly,  those  missions  with  hard  X-ray  imaging  capabilities 
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like  NuStar,10  Astro-H  (previously  NeXT),* 11  and  Simbol-X.12 
If  the  IC  emission  is  localized,  then  our  joint  fitting  technique 
can  be  used  for  many  much  smaller  regions  where  temperature 
mixing  will  be  less  significant  and  the  IC  component  will  be 
relatively  stronger.  Also,  the  B  field  strength  can  be  derived 
spatially  across  a  cluster,  yielding  a  better  estimate  of  the 
possible  pressure  support  provided  by  relativistic  components 
in  the  ICM,  which  could  modify  mass  estimates  that  depend  on 
the  hydrostatic  equilibrium  state  of  the  cluster  gas. 
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